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ANALYSIS OF SEVERAL RELATIONS 
AMONG ATMOSPHERIC STATISTICS 


By R. A. Minzner* and P. Morgenstern^ 

SUMMARY 

A statistically rigorous equation relating p, T , and p, the 
mean values of pressure, temperature, and density, respectively, 
into a gas-law-like expression includes (in addition to the regu- 
lar gas-law- like terms) a correction term proportional to the 
covariance between density and temperature. The application of 
the data from a set of 437 rocket soundings and several quasi- 
homogeneous subsets of the same data showed the value of the 
covariance term to be generally in the range of +1 to -1 per- 
cent of p for altitudes of 30 to 60 km. Covariance term values 
increased at greater altitudes and were largest for arctic 
latitudes and winter months. The application of the same data 
sets to another statistical equation relating 9a{p}/3z (the 
derivative of the standard deviation of pressure with respect to 
altitude) to a{p} (the standard deviation of density) in a form 
suggestive of the hydrostatic equation (which relates 8p/9z to p) 
yielded very poor agreement, undoubtedly owing to large uncer- 
tainties in the observed values of 8a{p}/3z. An integral form 
of the same equation, however, served to predict the observed 
data rather well. 


I. INTRODUCTION 

It has been customary in the generation of model atmos- 
pheres during the past decade (refs. 1-7) to use the following 
procedure : 

1. Determine a mean altitude profile for a single thermo- 
dynamic atmospheric variable from available sounding 
data , 

2. Introduce this altitude profile into the hypsometric 
equation and the gas law (or equation of state) , and 
compute the model altitude profiles of other thermo- 
dynamic properties. 

Thus, while the model atmosphere appears to be generally 
consistent with the mean values of temperature, pressure, and 
density deduced from sounding data, the model was derived as if 
there had been but a single sounding. It may be demonstrated 
that the mean value of density and the mean value of pressure 

*Electronics Research Center, Cambridge, Mass. 
tWalden Research Corporation, Cambridge, Mass. 



at a particular altitude calculated directly from the basic data 
of a number of soundings would not necessarily agree with the 
corresponding values of density and pressure as deduced by the 
above model-generation procedure. The implication of this situa- 
tion is that while the simple form of the gas law quite accurate- 
ly relates the thermodynamic properties of the atmosphere for a 
single sounding, it does not in general rigorously relate the 
corresponding mean values derived from a number of soundings. 

The gas law (or equation of state for an ideal gas) relating 
pressure, p, temperature, T, and density, p, can be expressed as 


P 



( 1 ) 


where M is the mean molecular weight of air and R is the univer- 
sal gas constant. The point values for p, p, or T may be ex- 
pressed as the sum of the corresponding mean values, p, p, or T, 
and of the departures, p', p', or T', respectively, as follows: 

P = P + P* (2) 

p = p + p' (3) 

T = T + T' (4) 

Substituting Eqs. (2), (3), and (4) into Eq. (1), and per- 

forming the indicated multiplications produces a second form of 
the gas law which, like Eq. (1) , is also applicable only to a 
single set of conditions, i.e., only to a single sounding at a 
given altitude. 


p + p' = ~ [p * T + p*T'+ p ' • T + p ' *T ' ] (5) 

Applying the averaging operator to both sides of Eq. (5) and 
noting that 


T ' = p ' = p 1 = 0 

yields the result 

P = ^ [p*T + p 7 *T ' ] (6) 
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where the quantity p'T' is commonly termed the covariance between 
p and T. Thus, the procedure for developing model atmospheres 
described above is strictly valid only if the covariance is equal 
to zero. 

Equation (6) can also be expressed in the alternate form 

P = ^ [p*T + a{p}»a{T}*r{p,T}] (7) 

where a{p} and a{T} are the standard deviation of density and 
temperature, respectively, and r{p,T} is the coefficient of 
correlation between them. Since a{p} and a{T} are almost always 
non-zero, this equation shows that the stated procedure for 
developing model atmospheres is strictly valid only if the cor- 
relation between p and T is zero. Both, Eqs. (6) and (7) , are 
statistical forms of the gas law as derived by Buell (ref. 8) . 

The covariance term in Eq. (6) represents the correction 
necessary to make the gas law strictly valid when applied to mean 
values for the thermodynamic properties. Stated another way, the 
covariance term is a measure of the error incurred during the 
development of model atmospheres by the method described earlier. 
Normalizing this error relative to the total pressure and con- 
verting the ratio to a percentage, the error may be expressed as 


percentage contribution of covariance term = 100 (8) 


Further insight into the factors contributing to this error term 
may be obtained by rearranging terms in Eq. (7) to produce 

p = ! P-T » g l - " «r{p,T}j (9) 

Thus the error term is the product of the coefficient of correla- 
tion between density and temperature, and of the coefficients of 
variation* of these two properties. The sign of the correction 
term is determined by the sign of the correlation coefficient 
since the other two quantities are always positive . 


*The coefficient of variation is defined as the ratio of the 
standard deviation to the mean value. 
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Another model relating statistical properties of atmospheric 
thermodynamic variables uses the hydrostatic equation as its pri- 
mary basis. The standard form of the hydrostatic equation is 



-g*p 


(10) 


where g is the acceleration of gravity and z is geometric altitude. 
It follows from Eqs. (2) and (3) that the hydrostatic equation 
also may be expressed as: 


'(ID 

From the basic statistical definition of the standard deviation, 
we have 


[a {p}] 2 = p' *p' (12) 

Differentiating both sides of this equation with respect to z 
produces 


0{p} i|iEi = (13) 

Substituting Eq. (11) into the right side of (13) yields 

a{p} = • g • p * • p ' (14) 


where g may be considered constant for a given latitude and given 
height*. By means of the definition of the correlation coeffi- 
cient, this equation can be expressed finally in the form 

= -g*a{p}*r{p,p} (15) 


where a{p} is the standard deviation of density as before, and 
r {p,p} is the correlation coefficient between pressure and 
density. 

*These restrictions could be removed by using geopotential, H, 
rather than the geometric height, z. This was not done in 
the present study. 


This is the expression derived by Buell (ref. 9) for the vertical 
gradient of the standard deviation of pressure at a constant 
height. It is analogous to an expression derived by Stidd 
(ref. 10) for the vertical gradient of the standard deviation of 
height on a constant pressure surface. 

The validity of Eqs. (6) and (15) was examined for the 30 
to 200 km region of the atmosphere by Minzner and Morgens tern 
(ref. 11) . The data collected for this purpose consisted of 437 
upper-air soundings covering the period 1947 to early 1965. 

These basic data were collected from 48 different sources in- 
cluding journal articles, institutional reports, and private 
communications. The data sample includes a variety of measure- 
ment techniques (falling-sphere, search-light probe, rocket- 
grenade, pitot-static probe, etc.), and covers a world-wide dis- 
tribution of 25 launch sites. The data set is considered to be 
a unique collection of upper-air soundings; it has been published 
in a separate technical report by Minzner, Morgenstern, and 
Mello (ref. 12) . 

The 437 soundings were stratified into quasi-homogeneous 
groups by a diurnal, seasonal, and latitudinal classification 
scheme. The samples within each of these subsets were used to 
validate the two statistical models, and disclose any observed 
time and/or space variations in the results. 

The model for the gas law was tested by evaluating the per- 
centage contribution of the covariance-correction term to the 
total mean pressure (see Eq. (8)). This quantity was calculated 
at integral kilometer altitudes for each of the diurnal-seasonal- 
latitudinal subsets of the data. The results were presented as 
an extensive set of 110\ machine-plotted graphs (ref. 11). 

The model based on the hydrostatic equation was tested by 
evaluating the expression: 


3a{p> 

percentage difference = 100 — — — 


[-g*a{p}*r{p,p}] 

3a{p} 

3 z 


(16) 


The statistical quantities a{p}, a{p}, and r{p,p} were calculated 
at integral kilometer altitudes for each of the diurnal-seasonal- 
latitudinal data subsets. The value of 3a{p}/3z was calculated 
using Stirling's central-difference formula (ref. 13) for 
numerical differentiation: 


3a{p} 

3z 


1 

Az 


f {z 1 >-f {z_ 1 > 
2 


f{z 2 >-2f{z 1 }+2f{z_ 1 }-f{z_ 2 } 

12 



(17) 
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where Az is the altitude increment between successive points (in 
this case 1 km), and f{zj_) is the value of a{p} at i altitude in- 
crements removed from the central value. A minimum of three 
points (one point on either side of the central value) and a 
maximum of five points (two points on either side of the central 
value) were used to evaluate the derivative. The results of the 
evaluation of Eq. (16) also have been presented as a set of 110 
machine-plotted graphs (ref. 11). 


II. DISCUSSION OF RESULTS FROM THE GAS LAW MODEL 

Variability of the Covariance Term in the Gas Law Model 

The graphs presented by Minzner and Morgens tern (ref. 11) 
are accompanied by a preliminary discussion of the physical 
implications. A more comprehensive review of the calculations 
and graphical results described previously has been conducted 
subsequent to the original publication. The results of this 
analysis are summarized below in terms of observed significant 
time and space variations. 

1. Height Variations 

a. Between 30 and 60 km, the magnitude of the correction 
term is usually much less than 1 percent, and only 
rarely exceeds 0.1 percent (to about 0.5 percent). 

b. The magnitude of the correction term tends to be greater 
for altitudes between 60 and 115 km than for altitudes 
below 60 km, but its value is still usually small com- 
pared with 1 percent for most of the stratified data 
sets. The number of data cells for which the values of 
the correction term exceed 1 percent for short altitude 
intervals increases with decreasing sample size at the 
higher altitudes. To increase the sample size at these 
altitudes, the data in different latitude belts were 
combined into a single subset. Combining these data 
introduces greater heterogeneity into the sample, re- 
sulting in increased correction-term values which grow 
to as much as 2 or 3 percent especially above 90 km. 
(Only for those data cells which present grouping of 
subsets across latitude bands is the sample size as 
large as 9 or more above 90 km.*) 

c. All the data were grouped into a single data set to pro- 
vide a sample size of 9 or more extending above 115 km. 


*A minimum sample size of 9 soundings was considered necessary 
to provide significant results. 
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In this case the value of the correction term increased 
to 5 percent between 115 and 155 km, whereas above this 
region its value increased to as much as 25 percent at 
210 km. 

2. Latitudinal Variations 

a. For the altitude region 30 to' 60 km, the value of the 
correction term tends to increase with increasing lati- 
tude. For the basic homogeneous data cells, this in- 
crease varies from < 0.1 percent in the tropics to 

< 0.5 percent in the arctic. 

b. For data cells combined across seasons and diurnal 
periods, the increased heterogeneity resulting from these 
compressions causes the value of the correction term to 
increase from 0.1 percent or less in the altitude region 
of 30 to 60 km in the tropics to values as great as 

2 percent for the same altitude region in the arctic. 

c. In the altitude region of 60 to 100 km, the trend is not 
so distinct since the number of available data sets is 
considerably less. The variation seems to be from values 

< 0.7 percent in the tropics to < 3.8 percent in the mid- 
latitude and subarctic belts. No arctic values are 
available at these altitudes. 

3. Seasonal Variations 

a. In the high latitudes, a strong seasonal variability is 
seen whereby the value of the correction term tends to 
increase from < 0.1 percent in summer to < 1.4 percent 
in winter, in the altitude region of 30 to 80 km. 

b. No significant trend in the variability of the value of 
the correction term as a function of season has been 
noted in the low and middle latitudes , where the value 
tends to remain < 0.1 percent for the same altitude 
region. 

4. Diurnal Variations 

a. Because of a lack of sufficient numbers of comparable day 
and night data, no definite trend can be said to have 
been detected in the variability of the correction term 
as a function of diurnal period. Only in the subtropi- 
cal summer data could such a comparison be made. This 
single comparison suggested a slightly larger value of 
correction term for nighttime data than for daytime data. 

b. When using data cells for which there has been a group- 
ing across seasons, the value of the correction term for 
tropical, daytime data showed no significant difference 
from that for tropical nighttime data. The subtropical 
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data combined across seasons reflected the same diurnal 
differences previously discussed for the ungrouped 
summertime data for that latitude. No similar compari- 
sons of seasonally grouped day and night data were possi- 
ble for the other latitude belts because of insufficient 
data. 

c. For grouping across both season and latitude, the value 
of the correction term for altitudes below 90 km is 
slightly larger for the nighttime data than for the day- 
time data; i.e. , < 1.0 percent compared with < 0.5 
percent. 

The observed time and space variations summarized above are 
based on analysis of the gross features of the numerical results. 
Within smaller altitude regions, values for the correction term 
showed one or more random negative excursions which are large 
compared with the remainder of that portion of the graph. A 
particular case in point is that illustrated by the graph for the 
tropical, summer, diurnal-mean data cell (see Figure 1). Sharp 
negative excursions of the correction-term value to -0.4, -1.4, 
and -2.4 percent are seen to exist at 51 , 58, and 68 km, re- 
spectively, whereas in the remainder of the altitude region 
between 35 and 65 km the term has values of < 0.1 percent. These 
appear to be the result of spurious density values at isolated 
altitudes and are most evident in data cells of small sample 
size. Consequently, these excursions have been largely disre- 
garded in the search for general trends . 

Analysis of the Covariance Term 

For a better understanding of the factors influencing the 
variability of the covariance term, it is convenient to review 
the second term of the bracketed factor in Eq. (9) . Of the three 
factors in this expression, the ratio a{p}/p as well as that of 
ct{t}/T vary only slowly with altitude. These two ratios have 
values which are greater than zero and less than +1.0 for all 
cases considered. The value of r{p,T} by definition may never 
exceed the range of -1.0 to +1.0. The sign of the triple product 
is determined by the sign of r{p,T), and the entire term vanishes 
when r{p,T} passes through zero. 

A more detailed review of the computed, but unplotted values 
of a{p}/p and a{T},/T, shows that both quantities tend to increase 
with altitude. The first quantity ranges from as little as 0.05 
or less at 30 km to 0.5 or more at 200 Jem. The second ranges from 
as little as 0.02 or less at 30 Jem to 0.4 or more at 200 km. 

Such changes with altitude are readily explained by a combination 
of greater relative variability in their parameters and increased 
measurement uncertainty of both p and T with increasing altitude. 
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TROPICAL 
Figure 1.- 


ALTITUDE (km) 

SUMMER DIURNAL MEAN 

Percent contribution of the correction term to the mean pressure 
of the Buell gas law equation as a function of altitude 



At any particular altitude, the values of both ratios tend to be 
smaller for homogeneous data than for heterogeneous data; i.e., 
smaller for data associated with a single latitude band, a single 
season, and a single diurnal period than for data which have not 
been segregated according to these variables. Because of these 
variations in the values of the two ratios , we might expect to 
find the values of the covariance or correction term to increase 
rather uniformly with increasing altitude and with increasing 
heterogeneity of other conditions of observation, provided that 
the value of r{p,T} were constant. The fact that the value of 
the correction term does not increase uniformly with altitude is 
due to a number of factors : 

1. The value of a{p}/p does not increase smoothly with 
increasing altitude. Spurious values of density increase 
the value of a{p} to a much greater extent than they 
influence the value of p. Hence, particular altitudes 

or limited altitude regions occasionally have abnormally 
large values of a{p}/p. 

2. Anomalous measured values of p at one altitude produce 
associated spurious computed values of T_at another 
altitude. Thus, a large value for a{T}/T at a given 
altitude may be a consequence of a large value of a{p}/p 
at another altitude. 

3. In addition, since single observations of p and T at any 
altitude are related by the gas law, a spuriously large 
value of p will have an associated spuriously small value 
of T, at that altitude. The converse also applies. This 
situation enhances the generation of the negative value 
of r{p,T} at the same altitude, and this value of r{p,T} 
tends to override the value which would otherwise be 
associated with that altitude region. 

4. Investigations beyond the scope of this study suggest 
that various factors cause the value of r{p,T} to vary 
in a somewhat systematic manner with altitude, latitude, 
and possibly with season. It is sufficient to indicate 
at this point, however, that r{p,T} is far from being 
constant. 

The relative influence of the factors a{p}/p, a{T}/T, and 
r{p,T} on the values of the correction term depicted in Figure 1 
are illustrated in Figure 2. It is apparent that the triple pro- 
duct of the three factors and the implicit relationships are such 
that spurious local values of density produce strong negative 
values in the altitude profile of the value of the correction 
term. Assuming that such spurious density values are not truly 
characteristic of the atmosphere but merely represent sampling 
fluctuations, it appears reasonable to disregard such random ex- 
cursions in the graphs of the altitude dependence of the correc- 
tion term. 
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SAMPLE SIZE 



Figure 2.- Altitude profiles for the correction term 
value and each of the three factors which 
contribute to this value as determined for 
the tropical, summer, and diurnal-mean 
data cell 


Summary of the Analysis of the Gas-Law Model 

In summation it appears that the absolute magnitude of the 
correction term of the statistical form of the gas law increases 
with increasing altitude and with increasing latitude, but gener- 
ally it remains below 1 percent for altitudes below 80 to 90 km, 
when a significant sample size exists. The value of the correc- 
tion term seems to be greater in the winter than in the summer 
for high-latitude data, but no such variation is seen for low- 
latitude data. Time of day does not appear to influence strongly 
the value of the correction term in our limited data sample. 
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Over broad height intervals below 60 or 70 km the correction- 
term values for homogeneous data remain near 0.1 percent. This 
value is similar to that reported by Wood and Spreen (ref. 14) 
for homogeneous data at altitudes below 30 km. 

For most model-atmosphere requirements at altitudes below 
120 km, the value of the correction term which should be 
applied in computing mean pressure is usually small compared with 
other uncertainties in the models, and the application of this 
correction term may be omitted. Consequently, for altitudes 
below 120 km, the procedure commonly used to develop model 
atmospheres appears to be a valid one. Sophisticated models 
based on observed data, however, might well be based on calcu- 
lations which include such a correction term, particularly for 
altitudes above 120 km. 


III. DISCUSSION OF RESULTS FOR THE STATISTICAL FORM 
OF THE HYDROSTATIC EQUATION 

Description of the Graphical Data 

The validity of the statistical model (Eq. 15) as applied to 
numerical data was determined by evaluating Eq. (16) . This 
quantity was found to be highly variable with changes of one or 
two orders of magnitude between adjacent values accompanied by 
frequent change of sign. The unconnected dots in Figure 3 illus- 
trate this variability for a typical data cell. This function is 
associated with the left ordinate of the graph which is designated 
PERCENT DIFFERENCE. Any computed values beyond the indicated 
range of the ordinate, ±100 percent, have been truncated at those 
limits in the plotting procedure. 

To improve the readability of this and similar graphs, only 
an envelope of the successive maxima and minima of the percent 
difference was plotted as a function of altitude. The graph of 
such an envelope for the same data cell depicted in Figure 3 is 
represented by the pair of dotted curves in Figure 4 . 

The series of + signs in each of Figures 3 and 4 represent 
the graph of a{p} as a function of altitude. The ordinate for 
this graph is given at the right-hand side of the figure. The 
scale is logarithmic over the range 10^ to 10"^ nt/m^ and the 
general linearity of the curve shows the near exponential height 
dependence of the statistic. 
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TROPICAL AUTUMN DIURNAL MEAN 


Figure 3.- Altitude profiles of the standard deviation of 
pressure and of the percent difference defined 
by Eq. (16) for the indicated latitude belt, 
season, •’nd diurnal period 
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Figure 4.- Altitude profiles of the standard deviation of 
pressure and of the envelope of the extreme 
values for percent difference using the same 
data set as in Figure 3 
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Interpretation of the Percent Difference Validity Test 

A review of the results given in Figure 4 and in the many 
similar graphs depicting the same model for other data sets, as 
presented by Minzner and Morgenstern (ref. 11) suggests that the 
agreement between the two sides of the statistical model is 
rather poor. The truncation procedure used in the preparation 
of these graphs is such that in many instances the true differ- 
ence is greater than the graphs suggest. However, it should be 
noted that the greatest departures occur in the regions where 
the curve of a{p} shows rapid changes of slope, resulting in poor 
agreement with the model. In height regions where the profile of 
a{p} is nearly linear on this graph, i.e., the derivative is 
nearly constant such as the 45 to 60 km interval, the agreement 
is much more favorable. 

A better indication of the relative smoothness of the indi- 
vidual terms in this model is shown by the curves in Figure 5. 

The crosses represent the values of the derivative on a logari- 
thmic scale* while the solid line in Figure 5 represents the 
values of the right-hand side of Eq. (15) . A visual comparison 
of these curves shows the greater variability in the derivative 
term. However, both curves are reasonably smooth and in relative 
agreement over that portion of the altitude range for which the 
sample size exceeds 20. 

The influence of greater sample size is illustrated by the 
results shown in Figure 6 for a more populous data cell. The 
region of reasonably good agreement extends over a much larger 
altitude interval in this instance than in Figure 5. Some rela- 
tively large differences are observed in Figure 6, primarily in 
the region of rapid decrease of sample size, and one positive 
value of the derivative is found in the same region. 

Figures 5 and 6 both show that neither side of Eq. (15) 
represents a very smooth function of altitude. The derivative 
term is highly sensitive to local changes in slope of the c{p} 
profile produced by measurement uncertainties in the data. This 
suggests that further smoothing of the data is desirable for 
evaluation of the differential-equation model in the form of 
Eq. (15) . Section IV contains a description of an alternate form 
of the model which provides implicit smoothing of the data. 


‘Because a{p} is a decreasing function with altitude, the deriva- 
tive generally is negative. Hence the absolute values have 
been plotted. 
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Figure 5.- Altitude profiles of the right and left sides 

of Eq. (15) evaluated for the tropical, autumn, 
diurnal-mean data cell 


Summary on the Degree of Validity 

The analysis of this model suggests that the number of 
soundings available for use in the current study generally is 
inadequate to establish the validity of the model. This is a 
consequence of the sensitivity of the statistical variables to 
sampling fluctuations in the data. The model yields the best 
results for those altitude regions containing the largest sample 
sizes . 
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IV. REVISION OF STATISTICAL HYDROSTATIC EQUATION MODEL 

The results of the previous section have shown that the 
statistical model of Eg. (15) is highly sensitive to random 
sampling fluctuations of the meteorological variables. The 
nearly opposite character of differentiation and integration 
suggests that an integral form of the hydrostatic-equation model 
may be developed to overcome the limitations of the differential 
form. Equation (15) can be integrated with respect to altitude 
to yield 

a{p} z = a{p} Q . - J z g*a{p} *r{p,p}dz (18) 


where the subscript ( Q ) refers to the uppermost altitude of the 
height interval of integration. The quantity a{p} 0 is the stand- 
ard deviation of pressure at height z Q . This model can be simpli- 
fied further if the acceleration of gravity is assumed to be con- 
stant. The validity of this assumption is examined in greater 
detail in Section V where it is established that the resulting 
error is not significant. 

Equation (18) with (g) assumed constant may be evaluated by 
use of an appropriate numerical integration technique. Because 
cr{p} is nearly an exponential function of height, a logarithmic 
integration formula is desirable. Minzner (ref. 15) has tested 
several simple integration techniques on atmospheric variables 
and found that the logarithmic trapezoidal rule yielded satis- 
factory results when applied to density profiles. Application of 
this rule to Eq. (18) produces 


a{p} ± 


a{ P } o 



(19) 


where 


Q = a{p} *r{p,p} 


and the subscript (i) refers to the number of integer kilometer 
increments below z Q . A more detailed discussion of this integra- 
tion method is given in Appendix A of this report. 

Equation (18) as approximated by Eq. (19) was evaluated 
numerically with the same set of data used to develop Figure 6. 
The observed value of a{p} at z Q is used as an estimate of 
a{p} Q . The effect of sampling errors in cr{p} 0 on the values 
of alp}^ as derived for the right-hand expression is dis- 
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cussed further in Appendix A. A comparison curve, the altitude 
profile of a{p} was calculated directly from observed values of 
atmospheric pressure. Figure 7 shows the results of these com- 
putations. These curves are significantly smoother than the 
results from the differential-equation model shown in Figure 6. 

A comparison of the two curves in Figure 7 shows that the dif- 
ferences are smaller than in Figure 6, and now tend to be sys- 
tematic in character. The discontinuities still exhibited by 
these curves generally are associated with sharp changes in the 
sample size of the data set. 

The validity of this model can be examined in a manner simi- 
lar to that used for the differential -equation model, namely: 


percentage 

difference 


100 


a{p> i - a{p> 

alpT 


( 20 ) 


(The direction of the difference in the numerator of Eq. (20) has 
been reversed to be consistent with the implicit signs in Eq. 

(16) introduced by negative slopes.) The altitude profile of 
Eq. (20) for the same set of data used in Figures 6 and 7 is 
shown by the dashed curve in Figure 8. The percent difference 
in this graph shows a maximum of 60 percent error with large 
segments of the profile indicating less than a 15 percent 
difference. A direct comparison of these results with those for 
the differential equation model also is provided by the solid 
curve in Figure 8. The departures in this latter case are seen 
to be substantially greater. It is apparent that the integral 
version of the model provides a more suitable basis than the 
differential model for investigating the existence of latitudinal, 
seasonal, or diurnal variations in the validity of the statistical 
model . 


V. CONCLUSIONS 

The objective of this study has been to evaluate the 
validity of two models for relating the statistics of atmospheric 
pressure, temperature, and density. These models were applied to 
data assembled from 437 upper-air soundings for the region 30 to 
200 km. The results of these calculations showed that the 
covariance term of Eq. (7) was generally in the range ±1.0 per- 
cent of the mean pressure for the altitudes 30 to 60 km. In this 
height interval the covariance was largest at upper altitudes, 
in arctic latitudes, and during winter months. Above 60 km the 
covariance term exceeded ±1 percent in many cases but these 
results are less reliable because of the more limited number of 
sounding available for analysis, particularly at the higher 
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jure 6.- Altitude profiles of the right Figure 
and left sides of the differ- 
ential model of Eq. (15) 
evaluated for the subtropical, 
autumn, diurnal-mean data cell 


STANDARD DEVIATION OF PRESSURE (nt/m* ) 



ALTITUDE (km) 


7.- Altitude profile of the standa: 
deviation of pressure compared 
with the same quantity inferred 
from the integral of Eq. (18) , 
using the same data employed i 
Figure 6 
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Figure 8.- Comparison between percent departure for the differential- 
equation model (solid line) and the model expressed by 
Equation (18) in integral form (dashed line) - corresponding 
values of standard deviation of pressure calculated directly 
from pressure data (crosses) 




latitudes. Because of a lack of sufficient comparable day and 
night data, no definite trend can be said to have been detected 
in the variability of the covariance term as a function of the 
diurnal period. 

The analysis of the model given by Eg. (15) suggests that 
the number of soundings available for use in the current study 
generally is inadequate to establish the validity of the model. 
This is a consequence of the sensitivity of the statistical 
variables to sampling fluctuations in the data. The model yields 
the best results for those altitude regions containing the largest 
sample sizes. Because of these limitations, a revision of the 
basic model has been proposed changing it from a differential 
equation into an integral form. The revised model was evaluated 
for a single case by means of a logarithmic trapezoidal integra- 
tion scheme. A comparison of these results with those from the 
original model suggests that the integral form provides a more 
suitable basis for future investigations of the existence of 
latitudinal, seasonal, or diurnal variations in the validity of 
the model. 

Gravity variations as a function of altitude and of latitude 
were considered in the initial processing of the individual 
soundings. In the subsequent analysis of the data, the results 
of many soundings were grouped according to geometric altitude 
both within a latitude band, and across latitude bands. This 
grouping at identical geometric altitudes across latitude bands 
is now seen to be a less desirable procedure, and has perhaps led 
to larger values of standard devitations of density and pressure 
than would have resulted if the data had been grouped according 
to equal values of geopotential (refs. 16, 17). 
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APPENDIX A 


NUMERICAL EVALUATION PROCEDURES 

The purpose of this appendix is to discuss qualitatively 
how the magnitude of the uncertainties in the model values for 
c{p) are influenced by the choices in each of two procedures 
in the evaluation of the expression: 


a{p} 



cr{ p} *r{p ,p} 


dZ 


(A-l ) 


The first of these procedures is the numerical evaluation of 
the integral term by means of a quadrature formula . Several 
quadrature formulas were considered from which the logarithmic 
trapazoidal rule was selected as providing a satisfactory 
approximation to the actual integral. The merits of some 
alternate choices are discussed later. The uncertainty which 
would be introduced into the value of a{p} by selecting a less 
suitable quadrature formula is small, however, compared with 
the uncertainty introduced into a{p} by the observational 
uncertainties in cr{p} *r{p,p} . 

An even greater uncertainty is introduced into a{p}, 
however, by a poor choice in the second procedural item, i.e., 
that dealing with the direction of integration , upward or down- 
ward. It will be shown that only downward integration yields 
acceptable values of percent uncertainty over the entire region 
of integration, although upward integration as implied by 
Eq. (A-l) may, in some instances, be satisfactory for a very 
limited region at the lower end of the range of integration. 

With upward integration the model values of a{p} are 
determined as a small difference between two large quantities, 
an intuitively undesirable procedure. It will be shown that in 
such a situation, the percent uncertainty propagated into the 
difference a{p} is an amplified weighted average of the percent 
uncertainties of both the terms comprising the difference. This 
propagated uncertainty in a{p} is never less than the percent 
uncertainty of the minuend, which is the integration constant, 
and grows without bounds as z increases.* 

With downward integration, however, the model values of 
a{p} are determined by a summation operation as in Eq. (A-2) , 
which comes directly from Eq. (A-l) : 


♦Both lower-case and upper case forms of Z used interchangeably in this appendix 
signify geometric height. 
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a{p}*r{p,p} dZ 


(A-2 ) 


a{p} = a{ p} Q 


/ 


In such an operation, it can be shown that the percent uncertainty 
in the sum, is a normal unamplified weighted average of the 
percent uncertainties of the addends, and consequently never 
exceeds the percent uncertainty of the larger addend. This will 
be seen to be the integral term for all but the first few kilo- 
meters of the range of integration. 

The weighting feature operates so as to minimize the effect 
of the percent uncertainties of small quantities involved in the 
sums or differences , even when these small quantities have very 
large uncertainties . This feature is most apparent and partic- 
ularly beneficial in the summation representing downward integra- 
tion, where the magnitude of the constant of integration a{p} Q 
is very small compared with that of the integral term over most 
of the range of an extended integration, but where the percent 
uncertainty of a{p} Q is large compared with that of the integral 
term. This large percent uncertainty in the integration constant 
is effectively eliminated from the uncertainty in model values 
of a{p} by the weighting feature. 

Both the weighting and amplification features which will 
be demonstrated analytically in the next section, favor the 
downward direction of integration. This conclusion is readily 
substantiated from an analysis of the error expressions associated 
with each of Eqs . (A-l) and (A-2). 


ERROR ANALYSIS 

To simplify the form of the expressions showing the error 
propagation into a{p} within each of Eqs. (A-l) and (A-2), it 
is desirable to simplify the symbolic notation of the terms of 
these equations prior to writing the error expressions. Thus, 
the quantity being computed, cr{p}, is changed to Q; the constant 
of integration cr{p} 0 is changed to C, and the entire integral 
term is designated by I. Then, using the Gaussian method we may 
relate the absolute uncertainties 6Q, 6C, and 61 in Q, C, and I 
respectively, through partial derivatives as follows: 


5Q . 6c) 2 + (|§ . 6i) 2 
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Because the variables C and I in Eqs . (A-l) and (A-2) appear 
only in the first power, in separate terms, and without 
coefficients, it follows that: 


6 Q = ^(6C) 2 + ' ( 6 1 ) 2 (A-4 ) 

Designating percent uncertainties in Q, C, and I as PUQ, PUC, 
and PUI , it follows that: 


6C = C 



1 

100 


C • PUC 


and 


61 = I 




PUI 


(A- 5) 


(A- 6) 


These forms are equally correct when applied to the percent 
uncertainties of both Eqs. (A-l) and (A-2). The conversion 
of 6Q to a form involving PUQ, however, differs in the two 
cases, and we have: 


<$Q = Q • ^ = 3^-q * (C-I) • PUQ (A-7 ) 


applicable to Eq. (A-l), and: 


<$Q = Q * = xJo * (C+I) * PUQ (A-8) 


applicable to Eq. (A-2). 

Combining Eqs. (A-4), (A-5) , (A-6) , and (A-7) into one equation, 

and then combining Eqs. (A-4), (A-5), (A-6), and (A-8) into 

another equation we obtain: 
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(A- 9) 


PUQ X = 


for upward integration, and: 


puq 2 = 




) 2 - 

(I 1 

• pui^ 2 

C 1 

Z 1 


and: 



V* 

(I 2 

• PUI 2 ) 2 

C 2 + 

X 2 



(A— 10) 


for downward integration. The subscript 1 added to the symbols 
Q, C, and I in Eq. (A-9) designates the particular values of 
these quantities applicable to Eq. (A-l) for upward integration. 

The subscript 2 employed similarly in Eq. (A-10) designates the 
particular values applicable to Eq. (A-2) for downward integration. 


NUMERICAL VALUES 

In order to investigate the error amplification as well as 
the weighting effects of Eqs . (A-9) and (A-10) with some degree 

of quantitativeness, it is necessary to have some reasonable 
estimates of the values of C^_ , C 2 / Iq and I 2 , as well as the 
uncertainties in and C 2 . These all follow from an estimate 
of the altitude profile of Q = a{p) and its uncertainty. The 
values associated with the lower and upper ends of such a profile 
serve respectively as reasonable values for the constants Cq and 
0-2 an< 3 their uncertainties. With these quantities we can derive 
estimates of the altitude profiles of Iq and I2 through the 
following relationships stemming directly from Eqs. (A-l) and 
(A-2) respectively: 


I 1 = C 1 - Q (A— 11) 


I 2 = Q - c 2 (A- 12 ) 


The number of upper-atmosphere observations available for 
analysis falls off rapidly with increasing altitude. The 
decrease occurs nearly stepwise in three height regions: (1) 

in the vicinity of 30 kilometers, (2) in the vicinity of 50 
kilometers, and (3) in the vicinity of 90 kilometers. At heights 
up to about 28 kilometers, the values of atmospheric pressure 
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and density, as well as their respective standard deviations, 
a{p} and a{p} are reasonably well known from many thousands of 
observations. For purposes of this discussion, the value of 
a{p} at an altitude of 28 kilometers, the low limit of the 
altitude profile is taken to be 1.45 x 10 2 nt m -2 , with an 
uncertainty of ± 6 percent. The same percent uncertainty is 
assumed for the value of a{p} which is taken to be 2.08 x 10 -3 
kg m“ 3 at this altitude. The 28-kilometer value of u{p} and its 
uncertainty have been adopted as Cj and SC^ respectively. 

Above 90 kilometers height, the number of atmospheric 
soundings decreases rapidly with increasing height. Between 
90 and 108 kilometers, the available data sample for the study 
reported in the main body of this paper for any particular 
height level, summed over all observation sites, decreased from 
84 to 37 samples. For this latter altitude, the value of a{p} 
was calculated to be 2.7 x 10 -3 nt m -2 . while the value of a{p} 
was calculated to be 3.69 x 10“-*- kg m -3 (ref. 11) . This value of 
a{p} has been adopted for the upper limit, of the altitude pro- 
file of a{p}, where the percent uncertainties of both a{p} and 
a{p} are assumed to be ± 100 percent. This value of a{p} and 
its uncertainty have been adopted for C 2 and 6C^ respectively. 

The 108-kilometer value of c{p} is seen to be five orders of 
magnitude smaller than the adopted 28-kilometer value of 
1.45 x 10 2 nt m“ 2 . It is apparent that a{p} or Q must vary 
approximately exponentially, i.e., in the manner of KC~^ Z , 
between these two altitudes. Using the cited end-point values 
of Q to define the function, it can be shown that the value of 
K is 6.504 x 10 3 nt m -3 , and that of A is 1.359 x 10 - ^ m - -*- . 

The value of a{p} is assumed to approximate a similar function 

ke -az such that the cited end points imply a value of 

9.592 x 10 -2 kg m -3 for k, and a value of 1.368 x 10~4 m“l for a. 


UPWARD INTEGRATION 


With the value of Cj and the altitude profile of Z 
established, Eq. (A-ll) yields the altitude profile of 1^, an 
approximation of which is shown in Figure A-l. For the sake of 
simplicity, no uncertainty band has been indicated for this 
quantity although the uncertainty in I is discussed below. The 
adopted value of is also shown in this figure, but in this 
case an uncertainty band of ± 6 percent is given. The value of 
Ci and its uncertainty band are given by three lines j, ^ 1 . 2 ' 
and Cj o, where Ci >2 represents the adopted value of 
1.45 x ‘lO 2 nt m“ 2 ,*and where- Cj j and C-. j represent, respec- 
tively, the + 6 percent and -6 percent departures from the 
adopted value. It is apparent from the figure that except for 
the very lowest portion of the range of integration, the value 
of Ij is of the same order of magnitude as the value of C-^. 
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For altitudes above the lowest 4 kilometers of the range of 
integration, the value of 1^ is greater than one half the value 
of Ci . Thus, for upward integration, it is apparent that except 
for the lowest altitude within the range of the integration, the 
magnitude of Cq and 1^ are such that their percent uncertainties 
both contribute significantly to PUQ, as expressed by Eq. (A-9) . 
For these altitudes, the relative magnitude of the two terms 
under the square-root sign depends primarily on the relative 
magnitudes of the percent uncertainties. 

At Z Q , the initial altitude, 1^ is zero and Eq. (A-9) 
reduces to: 


PUQ 1 = 




PUC 1 ) 


PUC. 


(A-13) 


As Z increases above Z , 1^ increases rapidly. At the same time, 
the growing value of (1-^ • PUI-^J^ increases the value of the 
numerator of Eq. (A-9) to some value greater than / (C^ • PUCq) ^ , 
while the growing value of 1^ decreases the value of the 
denominator to some value less than . The two effects occurring 
simultaneously cause the value of PUQ-^ to increase rapidly from 
its minimum value of PUCq. At an altitude of about 20 kilometers 
above the bottom of the range of integration, where 1^ approaches 
a value of about 0.9 of Cq, the value of PUQ^ becomes about 20 
times that of PUCn . As the integration proceeds to even greater 
altitudes above the initial level, the value of PUQ-^ approaches 
infinity as (Cq - 1^) approaches zero, essentially independent 
of the values of PUC^ and PUIq. It is apparent from Figure A-l 
that for a ± 6 percent uncertainty in Cq and no uncertainty in 
Iq, the value of (C^ - Iq) , and, hence that of Qq could become 
zero at an altitude Z x as low as about 53 kilometers. At 
greater altitudes, the values of Q would be meaningless. An 
uncertainty band of ± 6 percent associated with Ij would expand 
the point Z x to a region extending from a lower limit of about 
43 km to the top of the range of integration, 108 km. Within 
this region, the model value of Q-^ would be meaningless. 

From this analysis, it is apparent that upward integration 
yields reasonable values of Q-^ and PUQ-^ only at altitudes very 
close to the reference level. 
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Figure A-l.- Relative values of the two terms of Eq. (A-l) 
as a function of altitude for upward integra- 
tion from altitude z 


DOWNWARD INTEGRATION 

The introduction of the value of the integration constant 
C 2 i as well as the estimated altitude profile of Q into Eq . (A-12) 

yields an estimate of the altitude profile of I 2 r the integral 
term associated with downward integration. An approximation of 
this quantity is depicted in Figure A-2. As in the case of 1^, 
simplicity precludes the depicting of an uncertainty band for 
I 2 . For downward integration, from high to lower altitudes, 
the value of the integration constant C 2 is very small compared 
with the value of I 2 for much of the range of integration. The 
adopted value of C, = 2.7 x 10 -3 nt m -2 with an assumed uncer- 
tainty band of ± 100 percent, if drawn to scale in Figure A-2, 
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could not be distinguished from the abcissa of the graph. 

Though not apparent from the graph, it may be stated that the 
value of I 9 grows from zero to about ten times the value of C 2 
as the altitude decreases about 18 km from Z c . From this point, 
the value of I 2 continues to grow by a factor of about 10 for 
each successive 17— kilometer decrement in Z included within the 
integral. From this information, one may show that C 2 is very 
small compared with I 2 for all but the upper 20 to 30 kilometers 
of the range of integration. Thus, as the increasing range of 
altitudes encompassed by the integral becomes greater than 30 
km, the contribution of PUC 2 to PUQ^ through Eq. (A-10) becomes 
negligible even for an uncertainty m C 2 as large as 100 percent. 



Figure A-2.- Values of the integral term of Eq. (A-2) 
as a function of altitude for downward 
integration from altitude z q 


For this situation, (Z Q - Z) > 30 kilometers, Eq . (A-10) 

reduces essentially to: 


puq 2 = 


4 


d 2 • pui 2 ) 


= PUI. 


(A-14 ) 
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At altitude Z Q , where I 2 is zero, Eq. (A-10) reduces exactly to: 


puq 2 = 


{ 


(C 2 • PUC 2 ) 


= PUC. 


(A-15 ) 


so that PUQ 2 = PUC 2 , which is ± 100 percent for the adopted 
conditions. For altitudes within the upper 30 kilometers of 
integration, the uncertainty PUQo decreases from a maximum of 
PUC 2 at Z Q toward the value of PUI 2 as the altitude decreases 
from Z Q . Thus, Eq. (A-10) does not explode for any values of 
Z, and reasonable values of uncertainty of Q are obtained for 
the entire range of integration. 

Obviously, if only one direction of integration is to be 
selected, downward integration is to be preferred over upward 
integration. If the value of PUI 2 at 28 km altitude from down- 
ward integration is greater than PUCj, upward integration may 
possibly yield a slightly more precise value of Q in the immediate 
vicinity of 28 km. The uncertainty of such a value of Q will, 
however, rapidly grow to unacceptable values if this upward 
integration process is carried out for more than a few kilometers. 


QUADRATURE FORMULA 

A somewhat secondary but necessary consideration effecting 
the overall uncertainty of the model values of a{p} involves 
the selection of a quadrature formula to evaluate the integral 
term of Eqs . (A-l) and (A-2) . The selected quadrature formula 

serves as the basis for the development of an uncertainty 
expression for the integral term. The reason for designating 
the selection of a quadrature formula as a secondary consideration 
is due to the fact that while different normally acceptable 
quadrature formulae approximate the actual integral with varying 
degrees of accuracy, the error introduced by the poorest of 
these approximations is very small compared with other uncertain- 
ties. Three relatively simple formulas which have been considered 
are: (1) the linear trapezoidal rule, (2) the logarithmic 

trapezoidal rule, and (3) Simson's rule. Minzner (ref. 15) has 
shown the logarithmic trapezoidal rule to be a good approxima- 
tion of the definite integral for any function whose logarithm 
is nearly a linear function of the independent variable. Since 
the integrals 1^ and I 2 were both shown to meet this condition, 
the logarithmic trapezoidal rule has been adopted. Using this 
rule, we may rewrite Eq. (A-l) as well as Eq. (A-2) as: 
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g 


(A-16 ) 


a{p> i = cr{p} Q - 



where 

-2 

g = 9.80665 m sec , 
z = height in meters, 

V = cr{p}*r{p,p}, 

i = the number of increments included in the quadrature 

procedure thereby involving a total of (i + 1) altitudes. 


With downward integration, the value of (V^ - Vj_]_) is positive, 
the value of Zj_i is greater than Zj , sucnthat (Zj - Zj_^) is 
negative, and this negative sign coupled with the negative sign 
of the summation term effectively makes cr{p}j_ equal to the sum 
of two positive terms. The expression resulting from Simpson's 
rule would probably yield a slightly better approximation to 
the actual integral of Eq. (A-2) than that resulting from the 
logarithmic trapezoidal rule as shown in Eq. (A-16), but at 
the expense of considerably more complication. The linear 
trapezoidal rule would yield a somewhat poorer approximation, 
particularly in cases where the data points are distributed 
nonuniformly with altitude. The added uncertainty introduced 
into o{p} by the poorest of these three quadrature formulas 
would be small, however, compared with the uncertainty which 
is propagated into a{p} through the observational errors in the 
constant of integration and in the combined variable V, both 
of which enter into any quadrature formula that may be used to 
express I. 

The precise form of the equation used to express the 
uncertainty of the integral term depends directly upon the 
particular quadrature formula used, and while the form of the 
uncertainty expression related to each of the several quadrature 
formulae differs extensively, one from the other, the value of 
the computed uncertainty 6 I for any common set of conditions 
is not apt to differ greatly regardless of the form of the 
uncertainty equation from which it is computed. 

Again, the Gaussian method serves as the basis for the 
generation of the uncertainty expression, such that <51 the 
absolute uncertainty in the integral term is given by: 
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61 



For each related pair, Zj and Vj , of a data set, the un- 
certainty in both members of the pair may be lumped into one 
of the members Vj , with the other member Z . being considered 
exact. Consequently, the right-hand side of Eq. (A-17) may 
be simplified so that 61 may be expressed as: 



This expression implies the taking of the partial derivation 
of I with respect to each of the series of variables V Q , V-,, 

.... Vj_i, Vj , Vj +1 , •••» Vi . The quantity I is itself 

expressed as a series, formed in such a manner that each member 
of the series V Q to , except these two end values, appears in 
two successive terms of the series. Thus, in semi-expanded 
form, the integral term of Eq. (A-16) may be expressed as: 


1 = _g { (z i" z o) tln(v i /v o )] 1(v i“ v o ) + (z 2~ z i } t ln (v 2 /Vi ) ] 1 (v 2 -v 1 ) 
+ ... + ( Zj -z ) [ln(V j /V j _ 1 ) ]" 1 (V j -V j _ 1 ) 

+ (Z j+1 -Z j )[ 1 n(V j+1 /V j )]- 1 (V j+1 -V.) + ••• 


(A- 19) 

+ < z i-r z i-2> 

+ (Z . -Z . _!> [in (V^V^) ] - 1 (V i -V i _ 1 )} . 


The consequence of this expression is that the partial derivative 
of I with respect to V Q , and the partial derivative of I with 
respect to both result in expressions involving only one term 
of the series approximation of I . The partial derivative of I 
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with respect to all other members of the series of variables , 
i.e., V]_ to Vi_i results in e xpr essions involving two successive 
terms of the series approximation of I. Thus: 



The particular arrangement of the terms under the radical sign 
is to emphasize the similarity of two groups of terms. 

It has already been shown that, for the purposes of this 
study, a{p} may be assumed to be a simple exponential function 
of altitude varying from 2.083 x 10“3 kg m - ^ at 28 kilometers to 
3.687 x 10“8 kg m“3 at 108 kilometer altitude. Since r{p,p} 
usually has a value between 0.6 and 1.0, and never exceeds 1.0, 
we may approximate r{p,p} by unity of all altitudes, in which 
case : 


V = a{p} = ke az 


(A-21 ) 


where 

k = 9.592 x 10 
a = 1.368 x 10 


If the altitude increments of successive terms on the 
right-hand side of Eq. (A-20) are equal, the uncertainty 
expression can be considerably simplified since: 


V 


V. V.., V. 

2 _ _ 3 + 1 _ 1_ 

V. . V. " V. n 

3-1 3 l-l 


aAz _ . . . , , 

= e = a = 1.1466 


(A-22 ) 
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] 1 1 1 1 —QrAZ — — A O "7 0 0 

V. V. a - e = b - °* 8722 

D i 


In (V 1 /V Q ) = ln(V j /V j _ 1 ) = ln(V j+1 /V..) = ln(V ; ./V i _ 1 ) 


= In a = aAz = d = 0.1368 


With the above listed simplifications, Eq. (A-20) can be 
rewritten as: 



The radical on the right-hand side of this equation encompasses 
a series of (i + 1) terms, each of which contains an uncertainty 
factor 6Vj . This factor varies from term to term in accordance 
with term number and in accordance with the number of the members 
of the corresponding set of values of V-: comprising a part of the 
data set. In addition to the uncertainty factor, each of these 
terms includes a factor representing that portion of the altitude 
range to which the particular uncertainty 5Vj applies. For the 
first term of the series, corresponding to 6V 0 , the altitude 
factor consists of a single altitude element (Zy - Z Q ) multiplied 
by a dimensionless coefficient [ (a - 1 - d)/d2] which will 
hereafter be called F. The altitude factor of the last term of 
the series corresponding to 6Vy , also consists of a single 
altitude element multiplied by a different dimensionless coeffi- 
cient [ (b - 1 + d)/d2] which will hereafter be called L. The 
remaining terms under the radical , those terms corresponding to 
the general uncertainty factor SVj , all have altitude factors 
consisting of two terms. One of these terms consists of the 
altitude element (Zj - Zj_y) multiplied by the dimensionless 
coefficient L. The J other term of the altitude factor consists 
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of the altitude element (zj+i - z j ) multiplied by the coefficient 
F. These two altitude elements (zj - Zj_j) and ( z j+i “ z j ) are 
the two elements which are separated by Zj , the altitude 
corresponding to 6Vj . 

If (Zj - Zj_j_) is equal to (zj + i - z j ) , i.e., if all the 
altitude elements are equal to a fixed amount Az, Eq. (A- 25) 
reduces to: 


6 1 = g 



(L + F) 2 


3 =i+l 

E 


j=i 



2 




2 


(A-26 ) 


If the other assumptions, which lead to the series of Eqs . (A-21) 

to (A-24) , are applied to the coefficients defined as F and L, 
these coefficient may be defined in terms of a*Az: 


F = 



a • Az . n 

e - a • Az - 1 


(a* Az) 


(A-27 ) 


r = b - 1 + d 

J_, - ^ 


-a*Az , . i 

. e + a ♦ Az - 1 

(a • Az) 2 


(A- 2 8 ) 


Then, for AZ = 1000 m, a*AZ = 0.1368, a = e a ' Az = 1.1466, 
b = e -a* AZ = 0.8722, F = 0.5236, L = 0.478G, and F + L = 1.0016. 
If the values of the members of the set 6V 0 , •••, 6Vj , ••*, 6V± 

are known, 61 the uncertainty of the integral may be J calculated . 


For purposes of illustration, it is assumed that each member 
of the set. of 6Vj is 6 percent of the corresponding member of the 
set Vj . Thus, consecutive members of the set 6V 0 , ••*, 6Vj , ••*, 

6V^ form a geometric series with a common ratio "a" previously 
defined to be 1.1466. Eq. (A-26), however, calls for the sum 
of squares of the members of the set 6Vj , i.e.. 


j=i+l 

E 

j=i 
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The members of this set (6Vj) 2 also form a geometric series 
with the common ratio a2, such that: 


(F + L) 


T 

j=l 




(A-29) 


Evaluating this expression as a part of Eq. (A-26) on the 
basis of the specified conditions, yields a value 2.26 nt m - 2 
for 61, or 1.54 percent of 1.45 x 10^ nt m - ^, the 28 kilometer 
value of a{p}. It is interesting to note that for a sufficiently 
large number of increments in the integration, the percentage 
uncertainty of the entire integral is considerably less than the 
percentage uncertainty of any single element. This situation 
suggests that observations of Vj be made at reasonably close 
altitude intervals. 

It can also be shown that since the absolute uncertainty 
grows exponentially with decreasing altitude (for a fixed relative 
uncertainty) , the lowest 17-km region of the integration, within 
which the value of V changes by a factor of 10 contributes about 
99 percent of the total uncertainty in 61. Thus, even if the 
percent uncertainty in the values of Vj at the upper end of the 
range of integration in excess of 40 km should be found to be 
very large, these large percent uncertainties would contribute 
very little to the uncertainty in I at the lower altitudes. 


SUMMARY 

The following statements regarding computational procedures 
have been demonstrated either explicitely or implicitely in 
this appendix. 

(1) Only downward integration of VdZ , implying the sum 

of positive values of V and I is generally acceptable 
for computing a{p} over extended ranges of altitude. 

(2) Any of a number of reasonable quadrature formulas will 
introduce less uncertainty in a{p} than is introduced 
by the uncertainty in the observed values of V. 
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(3) The greater the number of observations of V within a 
particular altitude range, the smaller will be the 
uncertainty in the computed values of a{p}. 

(4) When the range of altitudes involved in a particular 
set of data of V vs Z , extends for 50 or more kilo- 
meters, the uncertainties in V at the high-altitude 
end of the range have little influence on the un- 
certainty in the computer values of a{p} at the low- 
altitude end of the range of observation. 
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